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Abstract. Time-resolved databases with large spatial coverage are quickly becoming 
a standard tool for all types of astronomical studies. We report preliminary results from 
our search for stellar flares in the 2MASS calibration fields. A sample of 4343 M 
dwarfs, spatially matched between the SDSS and the 2MASS calibration fields, each 
with hundreds to thousands of epochs in near infrared bandpasses, is analyzed using 
a modified Welch-Stetson index to characterize the variability. A Monte Carlo model 
was used to assess the noise of the variability index. We find significnat residuals above 
the noise with power-law slopes of -3.37 and -4.05 for our JH and HK S distributions 
respectively. This is evidence for flares being observed from M dwarfs in infrared 
photometry. 



1. Introduction 



Time domain studies of M dwarfs, using both photometry and spectroscopy, are a criti- 
cal tool for understanding the nature of magnetically driven flares.The frequency, inten- 
sity, and duration of flares are fundamentally tied to the star's ability to release magnetic 
energy generated by the internal dynamo. Flare events occur stochastically, and there- 
fore large samples are needed to shed light on the physical mechanisms which govern 
the flare rate, its evolution over time, and is variation with stellar mass. 

Seminal work was carried out by Lacy et al. ( |1976| ) in measuring the frequency of 
flares from M dwarfs which had been studied for several years in optical light (Moffett 
1974[ ). They found that flare occurrence rates, measured as frequency of flares versus 
energy of flare, could be described using a power-law with a slope between -0.4 and -1. 
Thus lower energy flares are seen more frequently. 

Recently, use of large area surveys has provided a fresh look at the rate, spatial 
distribution, and lifetimes of M dwarf flares. |Kowalski et al.| ( |2009| ) studied the rate 
of serendipitous flares found in repeat scans of the SDSS field known as Stripe 82, 
resulting in flare rates congruent with |Lacy et al.| ( 1976 1. Hilton et al. ( 2010) ) examined 
the repeat exposures from the SDSS DR7 spectral database, finding that flares are more 
commonly seen in the younger stellar population of the disk, and for later spectral types. 

Stellar flares exhibit a variety of lightcurve morphologies. Many last only seconds 
to minutes, with a characteristic fast rise and exponential decay, and relatively low am- 
plitudes. Others have been seen to evolve over many hours, with enormous increases in 
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Aux over a wide range in wavelength, and complex lightcurve evolution (e.g. Kowalski 
et al.|2010 



The next generation of large photometric surveys will be coming online in the near 
future. Work by Hilton et al. (these proceedings) predicts that M dwarf flares will be 
observed in nearly every frame of LSST observations, suggesting that flare stars will be 
a large source of contamination for other variability studies. 

Current photometric studies to measure rates of M dwarf flaring and magnetic 
activity have almost exclusively focused on bluer wavelengths. This is understandable 
as the flare continuum emission may be modeled with an 8,500K - 10,000K blackbody 
( [Hawley & Fisher 1992 Kowalski et al. 2010| ), and thus the greatest flux increases 
will be seen in the shorter wavelength optical bands. However, flares are known to 
produce emission over a broad wavelength range on the Sun, and it is possible that 
these events may also be observed in the infrared, albeit with reduced amplitude. Here 
we statistically characterize the stochastic increases in flux from M dwarfs at these 
longer wavelengths and attempt to identify flares. 



2. Sample Selection 

Our sample consists of spatially matched photometry from two different surveys. The 
2MASS Calibration Point Source Working Database (hereafter 2MASScal) consists of 
photometry from the repeat calibration scans of 40 tiles spread across the sky, tåken as 
part of the 2MASS survey (Skrutskie et al. 2006 ). This database is fully described in the 



2MASS Explanatory Supplement^] A calibration field was imaged multiple times every 
night to ensure repeatable and homogenous photometry for the main survey. 2MASScal 
provides between 562 and 3692 scans of each of the 35 main tiles. The five special tiles 
in the SMC and LMC were each scanned between 108 and 468 times. These databases 
are presently being mined, and follow-up study being conducted, for periodic variability 
from objects such as binaries (Davenport et al. in preparation). 

Previous work with the 2MASScal data has used broad JHK color cuts to select 
low-mass stars ( Plavchan et al.||2008[ ). To reliably determine the spectral type of low- 



mass stars, we cross matched 16 of the 2MASScal tiles to the Sloan Digital Sky Survey 
Data Release 7 (Abazajian et al.|2009 ) using a matching radius of 5". The left panel of 



Figure 1 shows the location of all the 2MASScal tiles, as well as the SDSS footprint. 
The 16 matched fields are shown as purple stars, yielding 32,044 individual objects 
with photometry in both SDSS and 2MASScal. 

The right panel of Figure 1 shows a color-color diagram for the SDSS-2MASScal 
matched data. For the 2MASScal filters, the mean magnitude for each object over all 
epochs is used. Note the large number of galaxies with red colors, a potentially ma- 
jor source of contamination for low-mass star studies. To ensure a clean sample of 
M dwarfs, we removed any point sources from this matched dataset which the SDSS 
OBJ TYPE flag did not call a star. Photometric spectral types were then determined with 
the covariance matrix technique presented in Kowals ki et al.| ( [2009| ) using SDSS r,i,z 
magnitudes. Because these stars are mostly nearby, the photometry was not corrected 
for interstellar extinction. Our final sample, having colors consistent with MO or later 
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Figure 1. (Left) The locations of 2MASScal tiles on the sky. A few notable fields 
have been labeled. The footprint of the SDSS is shown in grey. Fields with 2MASS- 
cal and SDSS photometry are shown as stars. (Right) Color-color diagram of the en- 
tire 2MASScal -SDSS matched sample consisting of 24,465 stars, and 7,283 galax- 
ies, as determined by the SDSS OBJ TYPE flag. The |Covey etaT1 ( |2007] l fiducial stellar 
locus is shown in black for reference. 



spectral type, contained 4343 stars, with ~ 3.8 x 10 6 2MASScal photometric measure- 
ments. 



3. Characterizing Variability 

Finding flares in such a large catalog of observations requires a robust determination of 
outliers in the stellar lightcurves. In this section we describe the quantities we used to 
find such epochs, and a model to estimate the reliability of our flare selection. 

3.1. The <D Statistic 

Many statistics are available to quantify the variability of stars. A popular one is the 



Welch-Stetson ind ex (|Welch & Stetson T993J). We adopted the modified Welch-Stetson 



index O used by Kowalski et al. (2009) which was designed to recover flares. Since 
2MASS observed the J, H, K s bands simultaneously, at a given epoch n, <1> jjj t „ is defined 
as: 

\ Cj,n j\ <T H ,n j 

where ra is the magnitude at the given epoch, ra is the mean quiescent magnitude, and <r 
the photometric error at this epoch. We also computed <&HK s ,n in the same fashion. The 
<!> statistic will be positive when both filters are simultaneously brighter or fainter than 
the mean value, and negative when one is brighter and the other fainter. We removed 
epochs where the star became fainter in at least one band to ensure that we did not 
recover eclipses or photometric errors such as clouds, but still provided negative O 
values. Specifically, we eliminated epochs with mjj > ra# for both distributions so that 
they would have the same number of epochs. 
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Figure 2 shows the $>jh and <$>hk, distributions. In both samples, as with the 
optical results of |Kowalski et al. ( 2009| ), the null distribution (O < 0) drops off for large 
values of O. The O statistic provides a measure of the flux ratio relative to the quiescent 
value, in terms of the errors for a given epoch. Flare epoch candidates are selected by 
defining a cutoff value of O that is characteristic of the noise. Epochs with O greater 
than this cutoff value, and an increase in brightness as described above, are probable 
flares. 




Figure 2. Histograms of the flare statistic <& JH (Left) and <S>hk s (Right). In both 
panels the black line is the data (<I> > 0) while the blue line is the null sample (<£ < 0). 
The model null <L> for each distribution, as described in §3.2, is shown in red. 



3.2. Modeling the Noise 

The null sample of the <I> distributions, shown in blue in Figure 2, should be represen- 
tative of the intrinsic statistical noise from the star from low-level stellar variability, 
or extrinsic sources. To model the noise, we created a Monte Carlo simulation of our 
dataset. For every star in the sample, random measurements were drawn from a gaus- 
sian distribution with a standard deviation equal to the mean photometric error for that 
object. This was done with the same number of epochs per object as our final catalog 
contained, thus reproducing the mean properties of the sample. 

The O statistic was then calculated in the same way as the data using Equation 1 
for both <&jh and ®hk s - Our resulting artificial O distributions were symmetric about 
cD = 0, as expected for uncorrelated noise between the two bands. We then made the 
same cuts as in our treatment of the real data to remove //-band flux decreases. The 
model null distributions had approximately half as many epochs as the model quiescent 
O > distribution due to these cuts. This too produced a symmetric O distribution, 
shown in red in Figure 2, thus allowing us to use the O < null sample as a realistic 
estimate for the quiescent ® > distribution. Any ® > measurements in excess of 
the null sample should be good flare candidates. 
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4. Results 



We subtracted the real null distributions, shown in Figure 2 in blue, from the data 
(<!> > 0). The cD residuals in Figure 3 show the number of epochs which have <1> values 
in excess of the well-characterized random noise. The residual distributions appear to 
follow a power-law distribution at high cD values, as shown in blue by the least-squares 
fit to the cD > 10 data in each panel. 

At O values less than about 10 the residuals turn over, deviating strongly from the 
single power law slope. This is the regime where the noise dominates, as shown in 
the model cp distributions from §3.2, and no complete census of very low amplitude 
deviations is possible. 




0.0 0.5 1.0 1.5 2.0 2.S).0 0.5 1.0 1.5 2.0 2.5 
log *jh log *HKs 



Figure 3. The residual O distributions (data - null) for the 2MASScal (Left) and 
LINEAR (Right) samples. The blue lines are rough power-law tits to the residual 
distributions. 



Our power-law fits have slopes of -3.37 ± 0.12 and -4.05 ± 0.07 for the ®y# 
and Q>hk s samples, respectively. We caution that these are rough estimates, as these 
residual distributions contain spectral types MO and later. It is known that the rate of 
flares and activity lifetimes vary across the M dwarf spectral sequence (West et aL 2008 1 
and ideally these fits should be done for each spectral type individually. The residual 
fits presented here are not the same as the power-law fits used in Lacy et al. (1976 1 
or |Kowalski et al.| (|2009[), which estimated the energy of the flare and the occurrence 



frequency. However, the trend of the residual distributions seen in Figure 3 suggest 
that the number of detectable flares decreases with increasing wavelength. This is not 
unexpected, since only the largest amplitude blue flares will have enough flux in the red 
and infrared bands to be detected. 



5. Future Work and Conclusions 

The large numbers of epochs for most of the sample (typically a couple thousand 
epochs per object) presents challenges in eliminating source confusion due to vary- 
ing observing conditions and seeing. Our preliminary investigations of objects with 
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nearest-neighbors closer than 2 arcseconds show that the 2MASS photometry in our 
database does not always properly separate objects. In the worst case, some "objects" 
show bimodal distributions in their (RA,Dec) positions. While rare, this can result in 
objects blending under poor observing conditions, and having abnormally high flux val- 
ues measured, which in turn can give false positives for flares, and may be the cause of 
the very low amplitude tail in the O residual plots. More stringent tests for the spatial 
matching and removal or de-blending of such objects are underway and will be pre- 
sented in future work with this dataset. Our future work will also include matching to 
the SDSS Stripe 82 time domain photometry, to provide variability analysis for stars in 
up to eight photometric bands. 

We have shown here that statistically significant transient luminosity increases in 
late type stars, with spectral types MO and later, are present in NIR survey photome- 
try. The power-law slopes become more negative at longer infrared wavelengths. This 
provides some of the first solid evidence for flares appearing in infrared bands, and 
suggests that flare detection is increasingly rare at longer wavelengths. 
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